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Abstract 

We study the lattice Nambu-Jona-Lasinio model with two flavors of Wilson fermions 
in the large N limit, where N is the number of 'colors'. For large values of the four- 
fermion coupling we find a phase in which both, flavor symmetry and parity, are spon- 
taneously broken. In accordance with general expectations there are three massless 
pions on the phase boundary, but only two of them remain massless inside the broken 
phase. This is analogous to earlier results obtained in lattice QCD, indicating that this 
behavior is a very general feature of the Wilson term. 
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1. INTRODUCTION 



Recently, Bitar and Vranas|l| presented an extensive study of the lattice Nambu-Jona-Lasinio 
model. This model, although interesting in its own right as an effective low-energy theory of 
the strong interactions, can serve as an ideal testing ground for the properties of the Wilson 
term. This is because the model can be easily studied in the large N limit, where in this case 
TV refers to the number of colors. In this letter we are concerned primarily with the nature of 
the symmetry breaking induced by the Wilson term. We will show that in the large N limit 
there exists a phase in which both parity and flavor symmetry are spontaneously broken. 
It is the existence of this phase transition that is responsible for the occurrence of massless 
pions in the model despite the fact that chiral symmetry is explicitly broken by the Wilson 
term. The same phenomenon is also responsible for the masslessness of the pion in lattice 
QCD||. The existence of this phase was missed by Bitar and Vranas. Nevertheless, most of 
our results in the symmetric phase are in agreement with theirs. 

The model under consideration is defined through the action 

S = EE(fW[%!/) + M\x,y)W(y) + ^[a 2 (x) + n 2 (x)}5(x, y)\, (1) 

x,y 4=1 I iV J 



where ip l (x) has also two flavor components in addition to the indices shown. The combina- 

N 



tion jt will be refered to as f3 below. The matrix M(x, y) is defined by 



M(x,y) = W 



(7 M - r)6(x + n,y)- (7u + r)6(x - fi, y) 



+S(x,y) 



4r + m + a(x) + ry 5 7f (a;) • f 



(2) 



In the above action auxiliary fields 7r and a have been introduced to decouple the four fermion 
interaction 10]. The terms proportional to the parameter r come from the Wilson term. In 
the following we will set r = 1. The Wilson term is of order 0(a) (a is the lattice spacing) 
in the naive continuum limit but nevertheless it has a pronounced effect on the theory. It 
gives the fermionic doublers masses on the order of the momentum space cutoff, O(-), and 
in QCD it also serves to produce the correct anomaly of the flavor singlet axial current in 
the continuum limit [[J. Note that the Wilson term explicitly breaks chiral symmetry. 

As it is well known, the Nambu-Jona-Lasinio model in the absence of the Wilson term 
breaks chiral symmetry spontaneously at a critical value j3 c . Above (3 C , the theory is chirally 
symmetric for m — 0, and < ^ip 1 > serves as an order parameter for the transition. With 
the Wilson term in place however, < ip 1 ^ 1 > is always non-zero and chiral symmetry is 
explicitly broken. Massless pions, if they exist in the model, therefore cannot considered to 
be the Goldstone bosons of broken chiral symmetry. We will show below that massless pions 
do exist and that they must considered to be Goldstone bosons of spontaneously broken 
flavor symmetry. In addition to the Goldstone pions there is one remaining pion which is 
massless only on the phase boundary and is associated with the spontaneous breaking of a 
discrete space-time symmetry, parity. In general, when there are rif flavors, there are nj 2 — 2 
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Goldstone bosons and one mode which is massless only on the phase boundary. The flavor 
singlet pseudo-scalar meson is always massive. 

In Sec. 2 we derive the formulas for the condensate and the masses in the large N limit. In 
Sec. 3 we present the results of our calculation, and finally Sec. 4 contains a short discussion 
of our results. 



2. LARGE N APPROXIMATION 



We start by integrating out the fermion fields in Eq. (1) obtaining the following effective 
action for the auxiliary fields: 

S = -— Tr(logM) + Tr(logM f ) + 2NPY, & 2 (x) + **{x) , (3) 
2 1 J x L 

where the trace extends over space, color, flavor and spin degrees of freedom. Since the 
action is proportional to iV we can proceed in the large N limit by evaluating the functional 
integral around the stationary point of the action. We write 

0(j(x] 07r(x) 

a(x) ~ a s + —+J-, t?(x) ~ 7r s e 3 + —+J-, N -> oo, (4) 



where we have accounted for fluctuations around the translationally invariant saddle point. 
Note, that we have allowed for the possibility that the isovector 7? might develop a non-zero 
condensate which we have arbitrarily chosen to point in the 3-direction in isospin space. 
Such a condensate of course breaks both, parity and flavor symmetry. 

In momentum space, the inverse propagator at the saddle point is given by 

M s {jp) = i 7m sinpf, + (4 - ^ cosp^) + m + a s + i^r^s- (5) 

The effective potential evaluated at the saddle point is 

Veff = -2Nn f J log[g(p)) + Nn f P(a\ + vr 2 s ) (6) 

where the function g(p) is defined as 

g(p) = J2 sin2 p^ + * 2 s + Hp) + m qf- ( 7 ) 

Following Bitar and Vranas[|IJ we define the constituent quark mass m q = m + a s and the 
function w(p) = 4 — J2fi C0S P^- Demanding that the linear term in the fluctuations around 
the saddle point vanish, leads to the two gap equations 

P f d A p a s + m + w(p) 

= ^2 -y — ^ — ' (8) 

s 2 J (2tt)^(p) 1 J 
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Obviously, n s = is always a solution to the last equation, but as it turns out it is not 
always the one that minimizes the effective potential. 

The inverse propagators are obtained from the quadratic fluctuations and are given by 

G-\k) = (10) 
P f d 4 p E M sin(p At + ^)sin(p M -^) + 7r s 2 -[ W (p + |) + m (? ][ W (p-|)+m (? ] 



and 



2 J (2tt)4 0(p + !)3(p_!) 



+ P f d 4 p E M sin(p At + %)sin(^-%)+7r s 2 + [^(p+|)+m g ][ W (p-|) + m (} 



2 J (2vr)4 g{p + h) g{p -k) 

In the last equation we have used that the pion propagator is diagonal in flavor space. Note, 
that G~}(0) = for a = 1, 2 due to Eq. (9) when tt s ^ 0. 

It is very difficult to calculate the mass, m CT , for the pseudo-scalar meson from the complex 
zeros of the inverse propagator in Eq. (10) [[L|. To obtain the mass of the pions, which is 
expected to be small close to the phase boundary, we can define the pion wave function 
renormalization constant Z n a and pion mass m n a such that 

hmG-l(k) = Z-}(k 2 + ml a ). (12) 



Unless Z^a becomes infinite, the pion masses can be computed from 

<* = ! im n ^L ( 13 ) 



Note that since parity is broken in the phase in which tt s ^ 0, the mixed ir — o propagator 
is nonvanishing. Physical states in this phase are not eigenstates of parity and are obtained 
by diagonalizing the mass matrix. 



3. RESULTS 

To make our results comparable with those presented in M, we choose N = 2 so that 
(3 = 2(3 in the following. We solved the gap equations (8-9) numerically on a 10 4 -lattice 
using a simple Newton procedure with the bare parameters ranging form < (3 < 2.5 and 
—2 < m < —6, each in steps of 0.01. We compared the results of the calculation where n s 
is allowed to take a non-zero solution to those results that are obtained for 7r s set to zero. 
We found throughout that n ^ usually minimizes the effective potential in Eq. (6), if such 
a solution exists. Some care must be taken in this calculation due to threefold solutions of 
the gap equation for a s when n s = and (3 < 0.75. In all cases we determined the solutions 
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for tt s and a s that minimize the effective potential. The results for a s and rc s from this 
calculations are plotted in Figs. 1 and 2. 

In Fig. 3, the continuous line represents the phase boundary between the regions where 
7r s = and the region where parity- flavor symmetry is broken and ir s ^ 0. Note that the 
region with tt s ^ disappears for (3 > 1.41. We have extended the calculation for m = —4 
far into the weak-coupling regime where (3 is large but there do not appear to be any more 
non-zero solutions for tc s . Along the dashed line in Fig. 3, the mass m q of the constituent 
quark vanishes. The intersection of the line where m q = with the phase boundary at 
(3 ~ 0.33 and m ~ —2.7 is the continuum chiral limit and corresponds reasonably well with 
the prediction of |lj]. 

Close to the phase boundary the pion masses are expected to be small. Thus, they can 
be calculated using Eq. (13). We find that the masses for all three pions are equal on the 
side of the phase boundary where 7r s = [as can be seen directly from Eq. (11)]. As one 
approaches the phase boundary we find that m^ 2 ~ m — m c . Although this is of course 
the behavior expected of a theory which breaks chiral symmetry spontaneously, the reader 
should understand that here it is simply a consequence of the fact that the critical exponent 
of the flavor-parity breaking transition has its mean field value at large N. In the broken 
phase, the masses of the ir 1 ' 2 remain zero, they are the Goldstone modes corresponding to 
the two unbroken generators of flavor symmetry. The mass of 7r 3 on the other hand is zero 
only at the critical point. It is to be considered as the inverse correlation length of an Ising 
like parity breaking transition. This behavior of the pion masses is summarized in Fig. 4. 
As we mentioned before, we were not able to obtain any useful results for the sigma mass. 
Its expected behavior is such that it is nonvanishing throughout the phase diagram||. 

Finally, we would like to mention that we were unable to find massless pions for values 
of (3 > 1.41. This seems to contradict what was found in [ ij (see in particular Fig. 8 of this 
reference). On the basis of our interpretation of the phase structure of the model one does 
not expect a massless particle in that region since there are no phase transitions there. 



4. CONCLUSION 

We have analyzed the phase structure of the Nambu-Jona-Lasinio model on the lattice with 
Wilson fermions in the large N limit. We found that in analogy to lattice QCD with Wilson 
fermions there exists a phase in which flavor and parity symmetry are spontaneously broken. 
As opposed to what is expected to happen in lattice QCD, the region of broken symmetry 
is a simply connected, compact region, symmetric around m = —4. There appears to be no 
interesting phase structure at weak coupling. Since a complicated phase structure has been 
found at weak coupling of the lattice Gross-Neveu model (an asymptotically free four-fermi 
model in 2 dimensions) with Wilson fermions ||], the absence of such an interesting phase 
structure at weak coupling may be related to the fact that the Nambu-Jona-Lasino model 
is non-asymptotic free. We feel that the most interesting aspect of our results is that one 
is able to obtain both parity and flavor violation in a vector-like theory in a well defined 
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approximation scheme. Due to the presence of the Wilson term, well known theorems [BJ 
about the absence of such breaking do not apply. We hope that the results presented here 
will help to make it easier to accept the parity-flavor breaking scenario for lattice QCD as 
well. Also, in some of the recent proposals for putting chiral fermions on the lattice, one is 
forced to work in a region where m < m c (the broken region in our language) [fj. A good 
understanding of the nature of this phase is therefore very important. 
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FIGURE CAPTIONS 



FIGURE 1: ir s as function of the bare parameters m and (3, obtained from the gap equations 
in (8-9). The value of ix s is non-zero only in a near-circular region (see Fig. 3) between < 
(5 < 1.41 and —5.5 < m < —2.5 that is symmetric with respect to the axis where m = —4. 
Note that for small (3, ir s appears to develop a singularity similar to exp{ — (A + m) 2 / f3} / \ff3 . 



FIGURE 2: a s as function of the bare parameters m and (3, obtained from the gap equations 
in (8-9). The value of a s is antisymmetric with respect to the axis where m = —4, and passes 
continuously through zero on this axis for all values of f3 > 0. 

FIGURE 3: Phase diagram for the Nambu-Jona-Lasinio model in Eq. (1) with spontaneously 
broken parity-flavor symmetry. ix s is non-zero only in the near-circular region that is sym- 
metric with respect to the dotted line where m = —4. Everywhere outside of that region 
7i s = 0, and the theory maintains the symmetry. Along the dashed line the quark mass m q 
vanishes. The intersection of the line where m q = with the phase boundary where = 
is the continuum chiral limit of the lattice theory. 

FIGURE 4: Pion masses near to the phase boundary for some generic (3 = 0.58. For 
m > m c = —2.804, all pion masses are equal and non-zero. The pion masses vanish for 
m = m c , and, while m^z 2 becomes finite, the masses for ir^ 1 ' 2 ^ remain zero in the phase 
where parity- flavor symmetry is broken and m < m c . 
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